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Abstract 

We deal with an inverse problem arising in corrosion detection. We 
-^ ■ prove a stability estimate for a nonlinear term on the inaccessible portion 

^-< \ of the boundary by electrostatic boundary measurements on the accessible 
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1 Introduction 
> 

r | In this paper we deal with the inverse problem of determining the nonlinear 

iy-s i term f(u) in the boundary value problem 
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(1.1) 



where Ti and T 2 are two open, disjoint portions of d£l and r£> = d£l\ (Ti U T 2 ) 
and / is a Lipschitz function such that /(0) = 0. 
S^ i Such a problem with a specific choice of the profile of /, has been introduced 

^ ' and discussed in recent years by M.Vogelius and others for the modeling of the 

electrochemical phenomenon of surface corrosion in metals [7], [13], [20]. In the 
above boundary value problem, f2 represents the metal specimen, Ti represents 
the corroded part of the boundary, which is not accessible to direct inspection, 
F 2 is the portion of the boundary which is accessible to measurements and Td 
is the remaining part of dQ, which is assumed to be grounded. 
The inverse problem thus consists in the determination of / when one pair of 
Cauchy data u\r 2 , f^|r 2 is known for one non-trivial solution u to (1.1). 
It must be noted that the models of nonlinearity which have been discussed in 
[13], [20] arc of the form 

f(u) = A(cxp(au) — cxp(— (1 — cx)u)) 
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and it turns out that are such that in the direct problem the existence and the 
uniqueness of the solution are not granted. For this reason, we found it necessary 
to require some additional a priori assumptions. First an energy bound of the 
measured electrostatic potential u 



[ |Vw(x)| 2 < E 2 



Next, an a priori bound of the Lipschitz continuity of /, namely 
\f(u) — f(v)\ ^ L\u — v\ , for every u, v £l. 



(1.2) 



(1.3) 



Moreover, we found it necessary to assume the knowledge of some additional 
information on the measured current density g on the accessible part of the 
boundary T2- More precisely, we assume a bound on the Holder continuity of 
g, with 



ll.9llc°-(r 2 ) ^ G 



(1.4) 



Also, we shall require a lower bound on the same current density g. Namely, we 
shall prescribe that, for a given inner portion r2 i 2r of T2 (see (2.10) below, for 
a precise definition) and a given number m > 0, we have 



HfflU~(r 2 ,2r ) > m > ° 



(1.5) 



We note that, indeed, a lower bound of this sort appears to be necessary. In 
fact, knowing the Cauchy data on T2 for a single potential u, one can expect to 
identify / only on the range of values taken by u on the inaccessible boundary 
Y\. Thus, as a preliminary step of the treatment of this inverse problem, it 
is necessary to evaluate the amplitude of the range of u on Ti, see Theorem 
2.1, in which we prove that the oscillation of u on Ti is bounded from below 
by exp(— (^) ), where c > 0,7 > 1 are constants depending on the a priori 
data only. Next, as the main result of this paper, we prove a stability estimate, 
Theorem 2.2, that is we show that if u±, u^ are two potentials corresponding to 
nonlincarities /i, f^, whose Cauchy data are close 

|]V>i -^2|U=(r 2 ) < e , 
ll.9i -52|U*(r 2 ) sje , 



where ipi 



ui L and — — 
112 dv 



gi with i = 1, 2 and if g = g\ satisfies the a 



priori bounds (1.4), (1.5), then the ranges of u\, U2 on Ti agree on a sufficiently 
large interval V such that 



length of V ~ exp 



On such an interval the nonlincarities /i, /2 agree up to an error ui(e) such that 



LO(S) 



log 



where < 8 < 1. 



The study of determination of nonlinear terms in elliptic inverse boundary value 
problems has produced various different results. Just to mention a few, Cannon 
[8], Bcretta and Vogelius [6], Isakov and Sylvester [12], Sun [18]. Among the 
specific peculiarities of the present study, in comparison with the mentioned 
previous results, we emphasize that the nonlinearity is part of a boundary con- 
dition and that existence and uniqueness for the corresponding direct problem 
are not assured, see, in this respect, Kavian and Vogelius [13] and Medville and 
Vogelius [16]. For the determination of a coefficient in a linear boundary con- 
dition, we refer for instance to Fasino and Inglcse [10], Alessandrini, Del Picro 
and Rondi [4] and Chaabanc, Fellah, Jaoua and Lcblond [9]. 

1.1 Plan of the paper 

In Section 2 we formulate our main hypotheses and state the main results, 
Theorem 2.1 and Theorem 2.2, whose proofs are deferred to Section 3 and 4 
respectively In Section 3 we prove some preliminary lemmas, in which we 
observe that, taking into account the a priori assumptions, the solution is Holder 
continuous with its first order derivatives in a suitable neighborhood of the 
inaccessible part r 1: (see Lemma 3.2, Theorem 3.3). In order to prove Theorem 
2.1 we need also a stability estimate near the boundary for a Cauchy problem 
(see Proposition 3.1). The proof of this Proposition is based on a method 
developed by Payne and Trytten [17], [19]. We conclude the Section with the 
proof of Theorem 2.1. In Section 4 we adapt the above mentioned stability 
estimate for the Cauchy problem (see Theorem 4.1), when both the Dirichlet 
and Neumann data are affected by errors in the L 2 norm, and at the same time, 
we obtain a stability result up to the boundary. In Proposition 4.2 we show that 
the solution is locally invertible along a suitable curve on Ti and we evaluate 
its length. Finally we prove Theorem 2.2. 

2 Main assumptions and results 
2.1 Preliminaries 

A priori information on the domain 

We shall assume throughout that f2 is a bounded, simply connected domain 
in M. n , n ^ 2 such that diam(Sl) ^ D and with Lipschitz boundary cWl with 
constants ro, M. More precisely, for every Xq E <9f2, exists a rigid transformation 
of coordinates under which, 

nnB ra {x ) = {(x',x n ) : x n > 7(0/)} (2.6) 

where x E R™, x = (a/, x n ), with x' E M rl ~\ x n E M and 

7 :£ ro (zo)cIR"- 1 ^M, 



satisfying 7(0) = and 



\h\\c°.i(B ro (x )) ^Mr , 



where we denote by 

n n I, I, \i(x)-i(y)\ 

hWc-HB^xo)) = \h\\L°°(B ro ( Xo )) + r- sup — — . 

x,y£B ro {z ) \* — y\ 

*¥^y 

Moreover, we assume that the portions of the boundary I\ are contained respec- 
tively into surfaces Si, i — 1,2 which arc C 1,a smooth with constants ro, M. 
More precisely, for a given a, < a < 1 and for any zq G Si, i = 1, 2, we have 
that up to a rigid change of coordinates, 

SinB ro (z ) = {(x',x n ) :x n = tpi{x')} , (2.7) 

where 

^ i :B ro (zo)cW l - 1 ^M 

are C 1,a functions satisfying y>j(0) = |Vy>,(0)| = and 

\\<Pi\\ci,<*(B ro (z )) < Mr , 
where we denote 

IMIci.«(B ro ( Zo )) = \\tph°°(B ro (z )) + r o||V<^|| L ^ (Bro(zo)) + 

+ r ^ sup : . (2.8) 

x,yeB ro (z ) \X-y\ a 

x=£y 

In particular it follows that if 

x G r.j and dist(.x , Tjj) > r , 
then 

Q,nB ro (zo) = {(x',X n ) € B ro (z ) : X n > fi(x')} , 

where ipi is the Lipschitz function whose graph locally represents dfl. Moreover, 
since J7 D B ro (zo) (~l Te> = 0, cpi must also be the C la function whose graph 
locally represents Si. We also suppose that the boundary of I\, within Si, is of 
C l,a class with constants tq,M, namely, for any zq G dTi, there exists a rigid 
transformation of coordinates under which 

dTi nB ro (z ) = {(a;', a; n ) G -S ro (2;o) : £ n = Vi(a;')) ^n-i = V'tfc")} 

where x' = (x",x n -i), with x" G M"~ 2 , x„_i G K and 

ipi :B ro (z ) CM"" 2 — >R 

satisfying -0i(O) = |V-0 4 (O)| = and \\ipi\\ci<<*(B T0 (z )) < M. 
We introduce some notation that we shall use in the sequel 

U [ p = {x G fi : dist(x, <9£7 \ r<) > p} , (2.9) 

r <iP = f/ p x n r< , (2.10) 



with i=l,2 and 

ftp = {x £ Q : dist(x, dn) > p} . (2.11) 

A priori information on the boundary data 

The current flux g is a prescribed function such that 

IMIc°.«(r 2) ^ G . 
A priori information on the nonlinear term, 

We assume that the function / belongs to C 0,1 (M,R) and, in particular, 

/(0) = and \f(u)-f(v)\ ^L\u-v\ for every u, v £ E . (2.12) 

We recall that a weak solution of problem (1.1) is a function u £ H l (fl), such 
that u\r D =0 and which satisfies 

f Vu-Vp= / gp+ [ f{u)p for all p € H 1 ^), p\ =0 (2.13) 

where u\ denote the trace on T^ . 

r D 

From now on we shall refer to the a priori data as to the set of quantities 
?*o, M, a, L, G, E, D, m. 

In the sequel we shall denote with rj(t) and u)(t), two positive increasing functions 
defined on (0, +oo), that satisfy 



7){t) ^ exp 



for every < t ^ G , (2.14) 



u(t) ^ C |log(f)H , for every < t < 1 , (2.15) 

where c > 0, C > 0, -f > 1, < 8 < 1 are constants depending on the a priori 
data only. 

2.2 Main theorems 

The statements of the main results are the following. 

Theorem 2.1 (Lower bound for the oscillation). Let fl,g satisfying the 
a priori assumptions. Let u be a weak solution of (1.1) satisfying the a priori 
bound (1.2) then 



oscu^n(\\g\\ L ^ {r22ro) ) 



where ?y satisfies (2.14). 



Theorem 2.2 (Stability for the nonlinear term /). Let m £ H 1 (0), i = 
1,2 be weak solutions of the problem (1.1), with f = /, and g — gi respectively 
and such that (1.2) holds for each m. Let us also assume that, for some positive 
number m, the following holds 

IISilU-(r 2 , 2ro ) > m > . (2.16) 



Moreover, let ipi = u, L , i = 1,2. There exists eo > onfo/ depending on the a 
priori data and on to smc/i t/iat, i/, /or some e, < £ < £o; we /lave 

||-0i -i'2\\mr 2 ) < £ > 
ll.9i -g2|U 2 (r 2 ) ^ £ , 



then 

where 

is such that 



||/i(m)-/2(u)|U~(vo <w(e) 
V = (a, /3) C [-OE, C*£] , 

p — a > 



2 

and r), co satisfy (2.14), (2.15) respectively. 

3 The lower bound for the oscillation 

Proposition 3.1 (Stability for a Cauchy problem). Let J7 satisfies the a 
priori assumptions and let v G H (O) be a solution of the following Cauchy 
problem 

Av = , in CI, 

v = if , on i\, (3 1?) 

dv 

— = h, onTi, 
av 

where ip,h G i 2 (ri) and the boundary conditions are considered in the weak 

sense. 

Then, for every Pi G Tx^r > v satisfies the following estimate 



\\ v \\L2(B p (p )nul n ) < C (\\<p\\ L 2 irip) + |]/i||i2(r 1>p ) +\\v\\m(n)) 



• (IMU 2 (r liP ) + ||^|U=(ri, p )) 

where p G { W il M i r o, A J[+ Ki i r o ) , Po = Pi + 4 vi+m^ r " -v,v is the outer unit 
normal to Cl at Pi and C > 0, < <5 < 1 are constants depending on p, r , n, M 
only. 

Proof. By a well-known estimate of stability for a Cauchy problem, we have 
that for any Pi G Fi, ro 

IMIz,2(,B p (Po)n[/i o ) < Ci (\\tp\\ L 2 {ri p) + ||Vw|| L 2 (rip) + |[u|| ff i(fj)) 

■ (IMU»(r ltP ) + ||V«|| L2(riip) ) <5 (3.18) 



whcrc p e (i7#iP ro ' i7& r °) and p ° = Pl + wrrm r ° ' "■ whcrc " is * hc 

outer unit normal to Q at Pi and C\ > 0, < rj < 1 arc constants depending 
on p, r , n, M only. For a proof we refer to [19], see also [17]. Let us define the 
function h G L 2 (dVt) as follows 

h(x), for a.e. x G I\ ro , 

1 



h(x) = < -TF~1 / ft ) for a ' C - x G r 2 . 

I 1 21 7ri, ro 
, otherwise . 



Let us consider the following Neumann problem 

f Az = 0, in O , 

— - = h , on oil . 
< ov 

Note that J 9n h = 0, hence a weak solution z G i7 1 (f2) exists and it is unique 
up to an additive constant. We select the solution z of (3.19) with zero average, 
it is well-known that for such a z the following holds 

IMIffi(n) < C 2 \\h\\man) ^ C3|WU*(r 1)ro ) 

where Ci and C3 are positive constants depending on the a priori data only. 
Let us set w = v — z, thus w solves the following Cauchy problem 

( Aw = , in a, 

w = ip- z , on ri iro , / 3 20 x 

dw 

— - = , on ri ro . 

By a standard boundary regularity estimate (see for instance [2, p.667]), we 

have that w e C l,l3 {U\ ) and the following holds 

2 r ° 

INIc^CC/ 1 , ) < C4||w||ffi(o) (3.21) 

where < (3 < 1 and C4 > are constants depending on r$ , M, a only. By an 
interpolation inequality, (see for instance [3, p. 777]) we have that 

2 r 

where C5 > and < 7 < 1 arc constants depending on M, a, ro only. 
Moreover, 



IM|j?i(n) ^ IMI-ffH") + lkllffi(n) < C 6 ^||u|| ff i(n) + ||^||^(r ltro )J (3.23) 
where C 6 = max{l,C 3 }. From (3.21), (3.22) and (3.23) it follows that 

\\\7w\\ L2{ri2ro) ^ C 7 (\\v\\ m{a) + \\h\\ LHrirn)) y~ (3.24) 

(iMI^OW + ||^IU3( ri , 2ro )) 7 • (3.25) 

Applying (3.18) to w and using (3.24) we obtain 

/ \ 1 — fq 

FlU=(B p (p )nc/ 2 1 ro ) < c (JMU 2 (rv ) + WHLHr^,) +\\v\\ H ^n)) 
■ (JMU'OW + IWU»(r liro )J • 
And the thesis follows with 8 = 777. □ 



Lemma 3.2 (Holder regularity at the boundary). Let u be a solution of 
(1.1), satisfying the a priori bound (1.2) then there exists a constant C > 0, 
depending on the a priori data only, such that 



l|w||i«>(Brp (a )nn) < CE, for every z e fi 



(3.26) 



\\u\\ c o, Hri) ^CE 

where < a < 1 is a constant depending on ro, M, n onfo/. 

Proof. From the weak formulation of the problem (2.13) we have 



(3.27) 



Vu- Vy> 



f(u)<fi 



where -Bio = Sio (zq) H f2 and Tio = dB^a. n Ti and y is any test function in 
i/ 1 (0) such that supp<^ C Bni U Tio . 
By (2.12) wc have that 



Vu- V<p 






and by a trace inequality (see [1, p. 114]) it follows that 



Vu-Vip 



* 



^ CL / \V(u<p)\. 



(3.28) 



From now on one can proceed by the standard iteration technique due to Moser 
(see for instance [11]) and (3.26) follows. By the local bound (3.26) and by ap- 
plying again the method by Moser leading to the Harnack inequality (see [11]) 
we obtain (3.27). □ 



Theorem 3.3 (C l,a regularity at the boundary). Let u be a solution of 
(1.1), satisfying the a priori bound (1.2), then u £ C l,a (G) and there exists a 
constant C p > 0, depending on the a priori data and on p only, such that the 
following estimate holds 



||«|| C i,- ( m) ^ C P E 



(3.29) 



where p € (0,rp) 



Proof. Since, by Lemma 3.2, we know that u £ C 0,Q (ri), by the Lipschitz 
regularity of / we have that 



du 
dv 



(x)=f(u(x)) £C»> a (T x ) 



By well-known regularity bounds for the Neumann problem (see for instance [2, 
p. 667]) it follows that u 6 C la (U p ) and the following estimate holds 



\u\\ c i, a(u i } C ||u||tfo,«(r ljf ) + 



dv 



oo.-fr,.) 



+ ||Vu|| L 2( n ) 



^ 



< C 



da 



dv 



E 



(3.30) 



C0.«(T lifi ) 



where C > depends on the a priori data and on p only. Moreover, we can 

estimate the C°' Q norm of — — in terms of E, in fact 

dv 



= sup 



du{x) 



<), 



p\ a 



sup 



c?ti(x) du{y) 

Bv Bv 



2/ ^ver, „ |^-y| Q 



9f 



L > 2 > 2 

By the Lipschitz bound (1.3) on / and by Lemma 3.2 we obtain 



. /n" !/(«(&)) -/(«(i0) I 

sup |/(«(a;))| + ( - ) sup 



On 



■ / m /P\" lu(#) — w(w)| 

^ L sup |u(z)|+.L(£) sup 



< L||u|| c o, Q(ri) < C£ . 
So inserting this estimate in (3.30) we have the thesis. 



\x-y\ a 



(3.31) 

D 



Corollary 3.4. Let u be as above, then, for every p > 0, the function S^ belongs 
to C 0,1 (Ti, p ), with Lipschitz constant L depending on the a a priori data and 
on p only. 

Proof. Let x and y be two points in T\ p then, by the assumption (1.3) and 
by Theorem 3.3, it follows that 



du(x) du(y) 



dv dv 



= |/(«(s)) - f(u(y))\ < L\u(x) - u(y)\ < 
< LC p E\x-y\. 



The thesis follows with L = LC p E. 



□ 



Proof of Theorem 2.1. Let e = osc u > 0, since u = on To, we have 
that 



ll«IU'(r ltro ) < Cie 



(3.32) 



where C\ is a positive constant depending on the a priori data only. By the a 
priori assumption (2.12) on /, we have that \f(u)\ ^ Ls, moreover, since 



du(x) 



dv 



\f{u(x))\ onf! 



then 

du 



dv 
By Proposition 3.1, it follows 



<|r liro |5L e . (3.33) 



\\u\\ L2{Bp{Po)nuLo) < C(e + Ef" s ■ e s (3.34) 

where C is a constant depending on the a priori data only. Since the boundary 
of Cl is of Lipschitz class, then it satisfies the cone property. More precisely, if 
Q is a point of dCl, then there exists a rigid transformation of coordinates under 
which we have Q = 0. Moreover, considering the finite cone 

f x ■ £, 

C = < x : \x\ < r , -j— j- > cos( 



with axis in the direction £ and width 28, where 8 = arctan^j, we have that 
C d Cl. Let us consider now a point Q € r2,r and let Qo be a point lying 
on the axis £ of the cone with vertex in Q = such that do = dist(QoiO) < 
^. Following Licberman [15], we introduce a regularized distance d from the 
boundary of Cl. We have that there exists d such that d £ C 2 (Ci) n C 0,1 (fj), 
satisfying the following properties 

dist(x,dCl) 

• 7o ^ tT\ ^ T 1 ' 

d(x) 

• |Vd(x)| ^ ci, for every x such that dist(x, dCl) < 6r , 

• ||d|| C i.« < c 2 r , 

where 7o,7i, Ci, c 2 , b are positive constants depending on M, a only, (see also [3, 

Lemma 5.2]). 

Let us define for every p > 

Clp = {x eCl : d(x) > p} . 

It follows that, there exists a, < a ^ 1, only depending on M, a such that for 
every p, < p $J aro, f2 p is connected with boundary of class C 1 and 

cip ^ dist(x, dCl) < c 2/ o for every x G 9fi p fl Cl (3.35) 

where ci, £2 are positive constants depending on M, a only. By (3.35) it follows 
that 

Using the notation introduced in the Proposition 3.1, we define the point P = 

P ° ~ aVi+m* r ° ■ v and Po = min ^ 32MVT+m r °' ^ sin61 ^ Moreover > let 7 be 

a path in CIpo. joining P to Qo and let us define {yi}, i = 0, . . . , s as follows 
"1 

2/o = Qo, Vi+i = l{U), where t t = max{i s.t. [-y(t) - 2/i| = 2p } if |P - y t \ > 2p Q 
otherwise let i = s and stop the process. 



10 



Now, we will use the three spheres inequality for harmonic functions (see for 
instance [14] or [5, Appendix E]) that is 



l-T 



u 2 < u 2 \ ■ 



B 3 PO (yo) \ J B PO {yo) J \ J B4 PO (yo) ) 

where < r < 1 is an absolute constant. Now since B po (yo) C B^ Po {y\) and 
since, by hypothesis ||u|| jyi (n) ^ E, then we have 



/ "'< / 

J B P0 (y ) \JB; 

An iterated application of the three spheres inequality leads to 



u 2 \ -E 1 
lB Po(yo) \ jB 3 P0 {yi) I 



f u 2 ^(f u 2 ) -E 1 

JB„ n (y„) Wb„„(:d.,) / 



Finally, since we have B po (y s ) C B 3M r (-PcOnt/^, then by the Proposition 

3.1 it follows 

/ u 2 ^C{(e + E) l - s -(e) s } TS ■ 

J B Po(vo) 

We shall construct a chain of balls B Pk {Qu) centered on the axis of the cone, 
pairwise tangent to each other and all contained in the cone 

„. f x-£ 

C = < X : \x\ < 7"o, -j — r- > cos( 



where 6' = arcsind 2 -). Let B po (Q ) be the first of them, the following are 
defined by induction in such a way 

Qk+i = Qk - (1 + v)pk£, , 
Pk+i = UPk , 
dk+i = pdk , 

with 

1 - sin 6' 

M ~ 1 + sin 0' ' 

Hence, with this choice, we have pk = fi Pa an( l B Pk+1 (Qk+\) C B^ Pk (Qk). 
Let us now consider the following estimate obtained by a repeated application 
of the three spheres inequality 

\W\\L*{B Ph {Q h )) ^ ll' i IU 2 (fl3p k _ 1 (Qfc-i)) < 

< H U llL 2 (B Pfc _ 1 (Q fc _ 1 ))ll w llL^(B 4pi _ 1 (Q fc _ 1 )) 

< C \\ u \\h(B m (Q )) < 

.C c{[(e + E) 1 - s -(e) s ] TS y k ■ (3-36) 



11 



For every r, < r < do, let k(r) be the smallest positive integer such that 
dk ^ r, then since dk = M "0> it follows 



|log(-n| IIob(t-) 

< k(r) ^ 



log^i 



logju 



1 



and by (3.36), wc have 

\H\mB Pkir)iQk W) ) < c{ [{e + E) 1 " 5 • (s) s f] 



fe(r-) 



(3.37) 



(3.38) 



Since, by hypothesis, T 2 is contained in a C 1,a surface and by the regularity 
assumption (1.4) on g, it follows, by the same argument used in Theorem 3.3, 
that u g C 1 ' (E/| ra ). 
Let x S r 2j 2r , a; G B p kM -i {Qk(r)-i), since w <E C l > a {Ul ') we have 



du(x) 
dv 


< 


du{x) 
dv 



C\x-x\ a C 



du(x) 



dv 



Integrating over B p k{r) -i {Qk(r)-i)i we deduce that 



du{x) 



dv 



^ 



^ 



^nl^-J JB p hir} _ 1 (Qkir)-l) 

2 

2 



du(x) 



dv 



'I -r 



: ( ^"1 



Wn(~V") JB PKr) x (Q k ( r )-i) 

Applying the Caccioppoli inequality, we have 



|Vu(:r)| 2 dx + 2C* 2 



4r 2 



du(x) 



dv 



< 



C 



(Pk-i) 



n+2 



'kW-iWw- 1 ) 



w(a;) 2 d.T + CV 



2a 



and since fc is the smallest integer such that d^ ^ r, then dfc_i > r, it follows 



du{x) 



di 



< ; . ^ A n+2 / U (x) 2 dx + Cr 2 

(rsmfl') -'Bp fc(r) _ 1 (Q fc (r)-i) 



From (3.38), we deduce 



du(x) 



dv 



<, 



i^{[(e + Ey- S -(s)rJ +Cr>° 



Let us define 



ct ( £ ) = [( £ + j e;) 1 - 5 .(e)*]", 

thus the previous inequality becomes 



du(x) 



dv 



C 



^(^)} T +Cr 



Now, using (3.37), we have 



do 



12 



where v = — log (-J logr. We have 



du (a 



ih 



^c 



a(e) 



Now minimizing the function on the right hand side, with respect to r, with 
r £ (0, ^), we deduce 



du(x) 



du 



CC[ lo 



1 



a(e) 



Since this estimate holds for every x £ I^^roi we infer 

du . „/, 1 

dv 



^ CI log 

i°°(r 2 ,2, ) V cr(e) 



where C is a constant depending on the a priori data only. Hence, solving for 
e, we can compute 

OK 



e ^ C exp < - 



dv 



L~>(T2,2 ro 



Note that, recalling the a priori bound (1.4), and choosing c = 2(1 — log CG 1 ) 
and 7 = ^t= one trivially obtains 



£ ^ exp 



, for every t g (0, G] 



a 



4 A stability result 

Theorem 4.1 (Stability for a Cauchy problem). Let Q, /j i = 1,2 and 

<?i satisfy the a priori assumptions described above. Let m £ iJ (f2), i = 1,2 
be weak solutions of the problem (1.1), wii/i / = ft and g = gt respectively and 
such that (1.2) holds for each Ui. 



Moreover, let tpi = m 



then, for every p G (0, ro) 



1,2. Suppose that 

||-0i -ih\\v>(Ti) <e , 
ll.9i - 92\\li(t 2 ) < £ > 



II wi -w 2 ||ai(r 1>p) ^w(e) 



(4.39) 



where to is given by (2.15) urai/i a constant C > which depends on the a priori 
data and on p only. 

Proof. Arguing as in Theorem 2.1, we find the following estimate 

IK -U2\\l"o(t. p) < w(e) , 

<w(e) 



9ui 9^2 
<9^ 3^ 



L °°( r i,«) 
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where a) is a positive increasing function of the type (2.15), such that 

ui(t) ^ C\ log(i)r 71 for every < t < 1 

where C > 0, < 71 < 1 are constants depending on the a priori data and on 
p only. By an interpolation inequality wc have 

||V t (ui - U2)||z,°o (ritp) < C\\ui - M2|li=o (r )||ui - w 2 ||ci°(r 1 . p ) 1 "^ 

- 2 

where /3 = -^-^ and C > depends on the a priori data and on p only, thus by 
Theorem 3.3 it follows that 

|[Vt(«! - « 2 )|U~(r ltP ) ^ G||«i - uall^^ sE 1 ^ 

where C > only depends on the a priori data and on p only. 

It follows that for every e < £q, with £0 depending only on the a priori data 

|[V(«l -U 2 )||L~( ri , p ) ^ 

dui du 2 
dv dv 



C 



+ ||V i (ui-u 2 )||L»(r llP )< 
i°°(ri, P ) 

< Cw(e)' 9 (4.40) 



where C > depends on the a priori data only. Hence, 

||ui - u 2 ||ic» ( r 1)P ) + ||V(m - «a)||L~(r llP ) < w(e) (4.41) 

where 



w(t) < Cw(£) p for every < t < 1 



□ 



Proposition 4.2 (Local monotonicity). Let u be a solution of (1.1) satis- 
fying (1.2), then there exist a point x £ I\ T and a direction £ S R™ -1 , |£| = 1 
such that, in the representation (2.7) ofY\ near x, the following holds 

\V x ,u{x\ <p(x')) • £| ^ 77 (||<?|U~(r 2)2ro )) , x' £ Us, = {x' = t-£ + x', \t\ < r} 

with 

r = min |^, ^^(|| ff || ioo(raaro) )| (4.42) 

where < a < l,di > are constants depending on the a priori data onfo/ anrf 
7] satisfies (2.14). 

Proof. Arguing as in Theorem 2.1, we can introduced a regularized distance, 
in the sense of Lieberman, on Si from the boundary of Ti and consequently 
construct connected sets I\ p for every p, < p sj &ro, which satisfy 

I\c 2 fc C fi ih C ri, Slh (4.43) 
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where < a < 1, C2 > C\ > are constants depending on M, a only. 
Since, by Lemma 3.2, u £ C 0,a (Ti), we have that by (4.43) it follows 



osc u Js _osc u ^ oscu — 2CE ( — ) c~2 Q 



Moreover by Theorem 2.1, we infer that 



_osc u^ V (\\g\\ L ^ (r2 ) )-2Ce(-^) 
r"i jl \ c i/ 



C2 



Possibly replacing c by a larger constant in (2.14) and taking 

ri =mm|? ? (||.g|| L o= (r22? . o) ),aciro,ro| 
we have that 

^osc u > ?KIMU~(r 2 , 2n) )) • (4-44) 

Let us set, for simplicity, r\ = rj ( ||<?||_L°°(r 2 2r ) ) ■ Since in the a priori assump- 
tions we have assumed that the portion Ti of the boundary is of C la class, 
then we can locally represent the restriction of u (the solution to (1.1)) to Ti, 
as a function of n — 1 variables, more precisely, for every xq 6 Ti, up to a rigid 
change of coordinates, we denote 

w{x') =u{x',ipi{x')) for all x G Y x f~l B ro (x ). (4.45) 

By (4.44), it follows that exist two points x and y in T 1 n, such that 

i] ^ w(x) - u(y) . (4.46) 

Let us consider a continuous path a C L x n joining x to j/ and let us define a 

sequence {xi}i = o,...,z as follows xq = x, Xi = a(si) where Si = max{i , \a(s) — 
Xi\ = ^} if \y — Xi\ > ^ otherwise let i = I and otherwise stop the process. 

(\ n— 1 
— ) , where C > is 

a constant depending on n only. 

Let us define 

Mj = m ax V t u(x) 
%(i,)nri 

4 

where V t denotes the tangential gradient on I\. Let M, i, x be such that 
x e Bn (x-) n Ti and 

M= max {Mi} = |V t u(x)|. (4.47) 

By (4.46) and the mean value Theorem, it follows that 

r] ^ \u(x) - m(xi) -\ h |w(xj) - u(y)\ < 

4 



< J2 M *T < ^ C i 



i=i, 
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where C\ > is a constant depending on the a priori data only. Thus we have 

M ^ -^- > . (4.48) 

Now we use the local representation of u as a function of n — 1 variables (4.45), 
within Ti n Bn (#?)• Let us define the direction £ = i v ^™i (x') . We shall further 
restrict the function w to the segment t ■ £ + x', with 

u(t) = w(i •£ + £') . 

Now, we look for a neighborhood Uq of t = such that 



K(*)l > 7^7" for ever y *e[/ . (4.49) 

It follows that for every \t\ < ^ 

|r/(0)-t/(t)|< 0,1*1° 

where C2 > is a constant depending on the a priori data only. 
Thus we have 

M=\v'(0)\^\v'(t)\+C 2 \t\ a . 



Hence by (4.48), 

Let us choose t in such a way 



^-C 2 \t\ a ^\v'{t)\ 



C2IC < 



2Ci 



Hence (4.49) holds with C/ = [-1", r], where r = min < ^-, ( ^^ ) >. The 



4 ' \ 2C1C2 

thesis follows, observing that v'(t) = gf = V X 'u(x',(p(x')) ■ £ and, possibly, 
by a further adjustment of the constant c in (2.14). □ 

Proof of Theorem 2.2. Let x G L iiTl , n, £ G R™ -1 be the point, the length 
and the direction introduced in Proposition 4.2, with u replaced with u±. Up 
to a change of coordinates, we assume £ = e\. Let 

Vi(t) = Uiit ■ £ + £, <pi(t- £ + #)) , i = l,2 , 

where x = (x' , tp\(x')) is the local representation of Pi near x. 
By Proposition 4.2 and assumption (2.16), we have that 

|«i(t)| *s T)(m) , for every t G U = [-n,ri] . (4.50) 

We shall denote by 7/1 = ?/(||<?i||z,°°(r, 2r ))■ By the stability estimate (4.39) of 
Theorem 4.1, we have that 

v' 2 (t) ^ 771 — w(e) , for every t £ Uo . 
16 



Thus choosing Sq such that 
we have 






Vi 



KWI > y i for every t E U 



(4.51) 



Thus the functions Vi are invertible on C/o, let us denote by Vi their respective 
images and by 



8 % :Vi->U , z = l,2, 



(4.52) 



their inverse functions. Let us observe that the intervals V\ and V 2 overlap on a 
sufficiently large interval V. In fact, by (4.50) and (4.51) it follows that Vi are 
monotone. Without loss of generality, let us assume they are both increasing. 
We have that, taken 

a = , o = — , 

2 2 

the following hold 

Vi(a) < Vi(t) < Vi(b) , for every t € (a, b) , i = 1, 2 . 
Moreover, since by the Theorem 4.1 we have 

IK -U2\\l°°(t, n) < W(e) 

' 2 

then, it follows that, for e < e , setting V = («i(a) + 2w{s),v\{b) — 2w(e)), for 
every u e V, there exists £ G (a, 6) such that v 2 (t) = it. 

Let us estimate from below the length of the interval V. By the mean value 
Theorem, (4.50) and (4.42), it follows that 

|«i(o)-t^(6)| = K(OI|6-o|>7iri . 

Thus the length C of V is bounded from below by 

C ^ Ti77i - w(e) . 

Hence, possibly adjusting the constant c in the definition (2.14) of r\, we have 
that 

C ^ 77(771) - lu{s ) > -77(771) > . 

Let us consider any value u £ V, then using the inverse function s % , we have 

u = vxis 1 ^)) = v 2 (s 2 (u)) . 
Let us estimate 



^-(s 1 (u)e 1 ,if 1 {s 1 (u)e 1 ))-—-(s 2 (u)e 1 ,ip 1 (s 2 (u)e 1 )) 



\h(u)-h(u)\ 
dv 
dv 
dv 



^-(s 1 (u)e 1 ,ip 1 (s 1 (u)e 1 )) - —-{s 2 {u)e 1 ,ip 1 {s 2 (u)e 1 )) 



^-(s 1 (u)e 1 ,p 1 {s 1 (u)e 1 )) - —-(s 2 (u)e 1 ,ip 1 {s 2 (u)e 1 )) 



dv 

du- A 
dv 
du 2 
~£h7 



< 
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where e x = (1, 0, • • • , 0) e R n_1 . By Theorem 4.1 it follows that, for all u e V, 

sC w(e) . (4.53) 



— (s 1 (u)ei,ipi(s 1 (u)ei)) - — (s 1 (u)ei,ipi(s 1 (u)ei)) 



< 



By Corollary 3.4, we infer that 

— (s 1 (u)ei, ipiis 1 (u)ei)) - -^j(s 2 (w)ei, v?i(.s 2 (u)ei)) 

£ (Is 1 ^) - s 2 (u)| + l^i^^ei) - ^i(s 2 (u)ei)|) ^ 
Z(l + Af)|s 1 (u) - s 2 (u)| . 

By the mean value Theorem, we find 

V2(s 2 (u)) = v 2 (s l (u)) + v' 2 (s)(s 2 (u) - s l (u)) 

where s is a point between s 2 (u) and s 1 (m). Since 

v 2 {s 2 (u)) = v 1 (s 1 (u)) , 

by (4.51) and by Theorem 4.1, it follows that 

\s\u)~s 2 (u)\ < l\ V2 ( s \ u ))- Vl (s\u))\^ 

m 

2 
^ — uj(e) , for every u E V . 

Finally, we infer that 

\fi{u) - f 2 {u)\ ^ u)(e) , for every u&V, 
possibly by a further adjustment of the constant C in (2.15). □ 
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